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random matrix theory for the epsilon regime are well reproduced. We conclude that improved 
staggered quarks near the continuum limit respond correctly to QCD topology. 
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1. Motivation 

The low energy regime of Quantum Chromodynamics (QCD) exhibits a rich and interesting 
phenomenology, including the C/a(1) axial anomaly, chiral symmetry breaking and the topological 
properties of the theory. There are a number of detailed predictions of the properties of these low- 
lying modes, such as the existence of an Index Theorem ^ ||] and the distribution of the first few 
eigenvalues in fixed topological charge sectors |g[ ^ § [7[ §]. 

Any correct discretization of QCD must agree with those predictions close enough to the 
continuum limit. Here we show that this is already the case for improved staggered quarks at the 
lattice spacings typically used in present-day simulations. In particular, we show that improved 
staggered fermions do respond correctly to the gluonic topological charge, and discuss why some 
confusion on this issue exists in the literature. Here we expand the work in [pi HIT]. For more 



detailed results see [11]. Related work has been presented in [12, 1 



2. Improved Staggered Dirac Operators 

The massless, gauge-invariant, One-link staggered Dirac operator on a d = 4 dimensional 
Euclidean lattice with spacing a is 

1 d 

Hx,y) = — £ [U^(x)d x+fi , y -H.c.] , 77 v (x) = (-l)Wx (2.1) 

with wo an optional tadpole-improvement factor given, in our case, by the fourth root of the mean 
plaquette. E) is antihermitian and obeys a remnant of the continuum 75 anticommutation relation, 
{B,e} = , with e(x) = ( — 1) 2 -^= 1 . As in the continuum case, its spectrum is therefore purely 
imaginary, with eigenvalues occurring in complex conjugate pairs, {=h'A s , X s G M} . The corre- 
sponding action describes N t = 4 "tastes" of fermions which interact via unphysical "taste break- 
ing" interactions that vanish in the continuum limit as a 2 . In such limit there is an SU (N t ) <S> SU (N t ) 
chiral symmetry, and the spectrum is therefore A^-fold degenerate. There is also an exact Index 
Theorem. At finite lattice spacing the chiral symmetry group is reduced to£/(l)<g)t/(l) and we do 
not expect to see this picture, for example there will not be an exact Index Theorem anymore. 

In addition to the One-link operator, we have also studied several improved staggered op- 
erators, the so called ASQTAD, Fat7xAsq and Hyp , designed to suppress the taste-changing 
interactions [0, [B], [T§]. 



3. Details of the Simulation 

We use a quenched, SU (3) gluonic action that is both tree level Symanzik and tadpole im- 



proved [J1V|, [18|, [19|, g0|, gJJ]. It differs from the one used by the MILC Collaboration |g2| only in 
small one loop radiative corrections. We have generated several ensembles of around 1000 config- 
urations. 

We measure the gluonic topological charge Q using two different cooling methods, to look for 
consistency. We discard the configurations for which both methods disagree. In the a = 0.093 fm 
ensemble fewer than 10% of the configurations have an ambiguous topological charge. For the 
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Figure 1: The positive half of a typical eigenmode spectrum for a configuration of Q = 2, in lattice units. 
The .x-axis is eigenvalue number. The Hyp operator gives results very similar to Fat7xAsq , and they are 
not shown for clarity. 



finest ensemble, a = 0.077 fm, this number goes down to around 2%. We stress that we only use 
cooling to determine the topological charge. All the Dirac spectrum measurements are done on the 
original thermalized ("hot") configurations. 



4. Index Theorem 

We begin our analysis by qualitatively comparing the low-lying modes of various staggered 
quark operators on a typical gauge background selected from those with Q = 2. In Fig. [I] we show 
both the eigenvalue and the chirality of the first sixteen (positive) eigenmodes. Near the continuum 
limit we expect to see first 2\Q\ = 4 near-zero modes with their chirality renormalised slightly away 
from unity. The remaining modes should have chirality near zero and divide into almost degenerate 
quartets. The spectrum looks quite continuum-like for all Dirac operators and we see a clear Index 
Theorem. For the improved operators we also see a very clear quadruple degeneracy in the non- 
zero modes. The renormalisation of the chirality away from 1 is small for the improved Dirac 
operators (around Z = 1.2), becoming as large as Z rs 4 for the One-link fermions. This result is 



fairly typical of operator renormalisation for staggered quarks [E3L 24 , 25 ] 



5. Comparison with Random Matrix Theory 

The chiral symmetries of staggered quarks are more complicated than for continuum QCD. 
At finite lattice spacing there are Nf = 16 pion states, and only one of these becomes massless 
in the chiral limit, the "Goldstone pion", with mass M G . The 15 remaining states have masses 
~ Mng arising from the taste breaking interactions. The (chirally extrapolated) masses are &(a ) 
and therefore zero only in the continuum limit. 

There are potentially two universal regions for different parameters of the system: 

£-regime: (Aqcd) -1 < L < (M NG y l , (5.1) 
e'-regime: (M NG y l < L < (M G y l . (5.2) 

In the first one the chiral symmetries are as in the continuum, corresponding to an effective theory 
of 16 massless pions. The low-lying non-zero modes have a near N r fold degeneracy, and follow the 
distributions corresponding to continuum QCD (and chiral lattice fermions [^, ^7, 28, 29, 30, 31]). 
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At finite lattice spacing there is a second, e'-regime, corresponding to an effective theory with 
a single massless pion. There is not even approximate restoration of the continuum symmetries, 
and the associated RMT has only U{\)<&U{\) chiral symmetry. The universal predictions for 
the eigenvalues are therefore strikingly different from those for continuum QCD. In particular, the 
predictions are the same for all sectors of topological charge p2[]. Presumably it is this regime that 



was studied in [J33J, [34J, |35|, pQ, [370. With the coarse lattices and unimproved gauge ensembles used 
in these studies, they observed no sensitivity of the eigenvalue spectra to Q, leading to the incorrect 
folklore that staggered quarks are "blind to the topology". To see the continuum chiral symmetries, 
there must be a sufficiently large mass gap between the heaviest pion and the lightest of the other 
hadrons. This requires the use of improved gauge and fermion actions and a sufficiently small 
lattice spacing. At coarse lattice spacings, with unimproved fermions, M^c Aqcd and there will 
be no £-regime. 

We separately fit the individual cumulative spectral densities to the predictions from RMT. 
These one parameter fits yield a prediction for the chiral condensate. We show some examples of 
the fits in Figs. |2[ It is clear that the improved Dirac operators match the predictions of RMT very 
closely over the full range of lattice spacings. By contrast, even at a = 0.077 fm the One-link 
Dirac operator still shows significant discrepancies. The fits also show clearly the expected varia- 
tion with topological charge. Each fit gives a value for the bare chiral condensate, and we find all 
of them to be broadly consistent with each other. 



6. Conclusions and Outlook 



We have studied the low-lying eigenmodes for the staggered lattice Dirac operator, using a 
range of improved and unimproved versions. The improved operators pass all the tests. In particu- 
lar, and contrary to previous accepted wisdom, such fermions do respond exactly as expected to the 
gluonic topological charge, and in a way identical to continuum QCD and other lattice formulations 
of the Dirac operator. 

We have seen that for improved operators the eigenvalue spectrum divides cleanly into near- 
zero and non-zero modes. The near-zero modes are characterised by a uniformly high chirality, 
with their relative number fixed by the Atiyah-Singer Index Theorem. The non-zero modes, by 
contrast, have chirality that is near zero and they divide into near degenerate quartets. In addition, 
the low-lying non-zero modes follow closely the universal distributions predicted by random matrix 
theory. 
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Figure 2: Comparison of the unfolded eigenvalue distribution with RMT for the lowest eigenvalue for 
different topological charges. The lattice size is kept at ah m 1 .5 fm 
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